Study of Xcj radiative decays into a vector meson 
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The decays XcJ — ^ {V ~ 4>, p", itj) are studied with a sample of radiative ip' 7XcJ events in 
a sample of (1.06±0.04) x 10* xj}' events coUected with the BESIII detector. The branching fractions 
are determined to be: ^ T-^) = (25.8±5.2±2.3) x 10~^ B(xci ^ 7P°) = (228±13±22) x 10"^ 

and B{xci — >■ 71^) = (69.7 ± 7.2 ± 6.6) x 10~®. The decay Xci — >■ 70 is observed for the first time. 
Upper limits at the 90% confidence level on the branching fractions for XcO and Xc2 decays into these 
final states are determined. In addition, the fractions of the transverse polarization component of the 
vector meson in xa — ^ decays are measured to be 0.29^0 {j-o'.og fo^' Xci — >■ "yp, 0.158±0.034to'.oi4 
for Xci 7P°5 and 0.247lg;Qg7lQ 026 ^'^^ Xci — ^ 70;, respectively. The first errors are statistical and 
the second ones are systematic. 

PACS numbers: 13.20.Gd, 14.40. Lb 



I. INTRODUCTION 

Doubly radiative decays of the type ip -fX — )• 77^, 
where V is either a 0, p'^, or lo meson, provide informa- 
tion on the flavor content of the C-even resonance X and 
on the gluon hadronization dynamics in the process 

. The spin and charge dependent couplings in radiative 
decays reveal detailed information which is particularly 
useful in the search for glueball and hybrid states Q ■ For 
the case where X = XcJ, the decay of the P-wave Xcj to 



"/V may provide an independent window for understand- 
ing possible glueball dynamics and validating theoretical 
techniques Q. 

Table U shows the theoretical predictions for Xcj radia- 
tive decays to a vector meson from perturbative quan- 
tum chromodynamics (pQCD) 0], nonrelativistic QCD 
(NRQCD) 01, and NRQCD plus QED contributions 
(NRQCD-hQED) [3|, and recent results from the CLEO 
experiment [6|. The experimental results for B{xci ~^ 
^p^ , 70;) are an order of magnitude higher than the cor- 
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responding theoretical predictions. However, by includ- 
ing non-perturbative QCD hadronic loop contributions, a 
recent pQCD calculation obtains results in agreement 
with the CLEO measurements. Improved measurements 
of Xcj radiative decays to vector mesons using the large 
BESIII ■)/'' sample will provide tighter constraints on the- 
oretical calculations. 

In this paper, we present measurements of radiative 
decays of the Xcj to the light vector mesons. The mea- 
surements have improved precision compared to CLEO's 
results, and Xci 70 decay is observed for the first 
time. In addition, the fraction of the transverse polariza- 
tion component of the vector meson in Xci "fV decay 
is studied, and the results indicate that the longitudinal 
component for Xci — ^ 7^ decay is dominant. This ob- 
servation may help in the theoretical understanding of 
Xci iV decays. 



TABLE I: Comparison of theoretical predictions on the 
branching fractions for Xcj radiative decays to a vector meson 
(in units of 10~^) and measurements from the CLEO experi- 
ment. The upper limits are at the 90% confidence level (C.L.). 



Mode 


CLEO [6] 


pQCD [2] NRQCD [3] 


NRQCD 










+QED [3] 


Xco -> ■yp" 


< 9.6 


1.2 


3.2 


2.0 


Xci 7P° 


243 ± 19 ± 22 


14 


41 


42 


Xc2 -5> 7P° 


< 50 


4.4 


13 


38 


XcO -i> 7^^ 


< 8.8 


0.13 


0.35 


0.22 


Xci — i> 'yijJ 


83 ± 15 ± 12 


1.6 


4.6 


4.7 


Xc2 — >■ 7^ 


< 7.0 


0.5 


1.5 


4.2 


XcO i4> 


< 6.4 


0.46 


1.3 


0.03 


Xci J<j} 


< 26 


3.6 


11 


11 


Xc2 7<^ 


< 13 


1.1 


3.3 


6.5 



II. BEPCII COLLIDER AND BESIII 
DETECTOR 

BEPCII/BESIII [1 is a major upgrade of the BESII 
experiment at the BEPC accelerator Q for studies of 
hadron spectroscopy and r-charm physics [l^l- The de- 
sign peak luminosity of the double-ring e^e~ collider, 
BEPCII, is 10^3 cm-^s-i at a beam current of 0.93 A. 
The BESIII detector with a geometrical acceptance of 
93% of 47r, consists of the following main components: 1) 
a small-celled, helium-based main draft chamber (MDC) 
with 43 layers. The average single wire resolution is 135 
/im, and the momentum resolution for 1 GcV/c charged 
particles in a 1 T magnetic field is 0.5%; 2) an elec- 
tromagnetic calorimeter (EMC) made of 6240 Csl (Tl) 
crystals arranged in a cylindrical shape (barrel) plus two 
end-caps. For 1.0 GeV photons, the energy resolution is 
2.5% in the barrel and 5% in the end-caps, and the po- 
sition resolution is 6 mm in the barrel and 9 mm in the 
end-caps; 3) a Time-Of-Flight system (TOP) for parti- 
cle identification composed of a barrel part made of two 



layers with 88 pieces of 5 cm thick, 2.4 m long plastic 
scintillators in each layer, and two end-caps with 96 fan- 
shaped, 5 cm thick, plastic scintillators in each end-cap. 
The time resolution is 80 ps in the barrel, and 110 ps 
in the end-caps, corresponding to a 2a K/tt separation 
for momenta up to about 1.0 GeV/c; 4) a muon cham- 
ber system (MUC) made of 1000 m^ of Resistive Plate 
Chambers (RPC) arranged in 9 layers in the barrel and 8 
layers in the end-caps and incorporated in the return iron 
of the superconducting magnet. The position resolution 
is about 2 cm. 

The optimization of the event selection and the esti- 
mation of physics backgrounds are performed through 
Monte Carlo simulations. The GEANT4-based simula- 
tion software BOOST 
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includes the geometric and 
material description of the BESIII detectors, the detector 
response and digitization models, as well as the tracking 
of the detector running conditions and performance. The 
production of the ip' resonance is simulated by the Monte 
Carlo event generator KKMC [l2[, while the decays are 
generated by EvtGen [l3| for known decay modes with 
branching ratios being set to the PDG [la] world av- 
erage values, and by Lundcharm (T6j for the remaining 
unknown decays. The analysis is performed in the frame- 
work of the BESIII Offline Software System (BOSS) ^ 
which takes care of the detector calibration, event recon- 
struction and data storage. 



III. DATA ANALYSIS 

A data sample of (1.06 ± 0.04) x 10* ip' events col- 
lected with the BESIII detector is used in this analysis, 
and an independent sample of about 42.6 pb~^ taken at 
^/s = 3.65 GcV is utilized to determine the potential 
background contribution from the continuum. In this 
paper, we focus on the exclusive decays of liXcJ, 
Xcj JhV, where 'jii'jh) designates the lower (higher) 



p , or w meson. The 



"■'^ and Lj candidates are reconstructed in the K~^K 



energy photon and V is either a 

TT+TT , and TT'^TT tt'^ decay modes, respectively. 

Charged tracks are reconstructed in the MDC, and the 
number of charged tracks is required to be two with net 
charge zero. For each track, the polar angle must satisfy 
I cos 6*1 < 0.93, and it must be within ±10 cm of the inter- 
action point in the beam direction and within ±1 cm of 
the beam line in the plane perpendicular to the beam. 
Since the efficiency of particle identification (PID) is 
lower for higher momentum (> 1 GeV/c) charged tracks, 
only the lower momentum charged track is required to 
be identified as a ii' or tt. 

Electromagnetic showers are reconstructed by clus- 
tering EMC crystal energies. The energy deposited in 
nearby TOP counters is included to improve the re- 
construction efficiency and energy resolution. Show- 
ers identified as photon candidates must satisfy fiducial 
and shower-quality requirements. The photon candidate 
showers reconstructed from the barrel region (| cos0\ < 
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0.8) must have a minimum energy of 25 MeV, while those 
in the end-caps (0.86 < | cos 9\ < 0.92) must have at least 
50 MeV. The showers in the angular range between the 
barrel and end-cap are poorly reconstructed and excluded 
from the analysis. To eliminate showers from charged 
particles, a photon must be separated by at least 10° 
from any charged track. EMC cluster timing require- 
ments are used to suppress electronic noise and energy 
deposits unrelated to the event. 

In order to choose the correct combination and improve 
the mass resolution, a four-constraint kinematic fit (4C- 
fit) is done under the assumption of energy-momentum 
conservation. Candidates with ^ 100 for this fit are 
retained. If an event has more than one candidate, the 
candidate with the smallest is kept. 

For ip' — > 77a; (oj — > Tr+TT^Tr") candidates, the gammas 
from the 7r° decay are selected as those that give the 
minimum of 



A/jr + TT 7172 



regions are defined as 1.05 < Mj^+i^- < 1.07GeV/c^, 
1.25 < M^+^- < 1.65GeV/c2, and (0.68 < M^+^-^o < 
0.71GeV/c2 and 0.85 < M^+^-^o < O.SSGcV/c^) for the 
(/), p", and w, respectively. 




IUl<7t3r) (OeV/c=) 
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'A^^l. 







where M^^^^ is the invariant mass of the photon pair, 
MttO (M^) is the nominal mass of tt" (oj), (Jt^q (^cj) is the 
mass resolution determined from MC simulation and is 
about 7 MeV/c^ ( 6 McV/c^). A tt" mass constraint for 
the 770; channel is included by doing a five-constraint 
kinematic fit (5C-fit), and events with < 100 are 
kept as JUhi^ (w — > Tr+Tr^Tr") candidates. 

To suppress background from multi-photon hadronic 
decays of the ip', lAf^,,^^ - > 25 MeV/c^ for 77^, 



> 600 McV/c2 for 77^", and 



MJ > 



25 MeV/c^ and |M^,^,^ - M^o| > 15 MeV/c^ for 77W are 
required. Here M^o and Af^ are the nominal masses of 
TT*^ and 1], respectively. The background from ip' — jhv' j 
rj' — > jiV (V = p'^, uj) is suppressed by requiring 
\M^,v-M^, \ > 15 MeV/c2. 
In XcJ 7/5*^1 there are potential backgrounds from 



QED e 



7e^e 



and /i^/i events where the leptons 



are misidentified as pions. To reject electrons, the ratio 
of the energy deposited in the EMC to the momentum 
measured in the MDC (£'emc/cpmdc) of tracks must be 
less than 0.8. To reject muons, tracks are removed if 
the number of layers with hits in the muon chamber is 
greater than three. The QCD backgrounds remaining can 
be effectively eliminated by requiring the opening angle 
between the two pions, cos^^+tt-, satisfy cos^tt+tt- > 
—0.8, and that between the two photons, cos ^7,7^ , satisfy 
—0.98 < cos 6*^,^^ < 0.5, in the laboratory frame. 

Figure [T] shows the K^K~, tt+tt", and Tr+Tr^Tr*^ in- 
variant mass distributions for the candidate events. The 
curves show the best fit to the mass spectra using 
a s-dependent Breit-Wigner function for signal and a 
polynomial for background. Events with \Mf^+j^- — 
M^\ < 0.01 GeV/c^, \M^+^- - Mp\ < 0.2GeV/c^, and 
|M^+^-^o-Af^| < 0.035GeV/c2 are taken as (ji, p°, and a; 
candidates, respectively. Here M^, Mp, and M^^ are the 
nominal masses of these vector mesons. The sideband 



FIG. 1: Invariant mass distributions of (a) K'^K~ , (b) tt^tt", 
and (c) ■K'^n~'K^. Dots with error bars are data; dashed lines 
are signal shapes; and dotted lines are the polynomial back- 
ground contributions. The signal regions and sideband re- 
gions are indicated with the solid and dashed arrows respec- 
tively. 

After applying the above criteria, there are still sev- 
eral peaking backgrounds in the XcJ mass regions from 
XcJ decays into non-7y modes with the same final states. 
From MC studies, the shapes of these backgrounds are 
found to be similar to those of the vector meson sideband 
background events. The invariant mass distributions of 
"fhV, where V ~ (j>, uj, respectively, are shown in 
Figs. [2Ka)-[2l^c). There are clear Xci signals in all decay 
modes, while Xco a-nd Xc2 signals are not evident. In or- 
der to extract the signal yields from the mass spectra, 
we first obtain signal shapes for each Xcj ^ jV mode (9 
decay modes in total) using MC simulations. Each of the 
distributions in Fig. [2] is fitted with a background shape 
composed of the vector meson mass sideband distribu- 
tion plus a 2nd order polynomial function and three Xcj 
resonances as the signal shapes. Parameters of the poly- 
nomial function and the normalization for each of the 
XcJ resonances are allowed to float in the fit. The fitted 
yields are summarized in Table IIIII Xci ^ 7P and 70; 
arc observed with a statistical significance larger than 
lOcr, and the significance for Xci — ^ 70 is 6Aa. Here, 
the significance is determined from ^— 21og(£o / ^mux) , 
where £max is the maximum likelihood value, and £0 is 
the likelihood for a fit with the signal contribution set to 
zero. Branching fractions are calculated after considering 
the signal efficiency, as listed in Table IIIII and the upper 
limits at the 90% C.L on the branching factions of Xrn 
and Xc2 decays are estimated by a Bayesian method [l8| . 
The effects of both the statistical and systematic uncer- 
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tainties to the upper limits are taken into account. All 
results are listed in Table IIIII 
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FIG. 2: Invariant mass distributions of (a) 'y(j), (b) -fp", and 
(c) 701. Dots witli error bars are data; histograms are the best 
fit; dashed histograms are signal shapes; and the gray-shaded 
histograms are tire sum of the sideband background and the 
background polynomial. 



IV. ESTIMATION OF SYSTEMATIC 
UNCERTAINTIES 

Tablc[n]shows the summary of all sources of systematic 
uncertainties. Many systematic uncertainties arc deter- 
mined using clean, high statistics control samples that 
allow results from MC simulation to be compared with 
those from data. 



A. Photon efficiency 

The photon detection efficiency and its uncertainty are 
studied by three different methods. 

The missing photon method uses a sample of ip' — > 
TT+TT" J/'0, J/i/j ^ /5°7r° events. Using events with four 
charged tracks, identified as pious, plus a good photon, 
the missing momentum is determined and used to predict 
the direction and energy of the missing photon. To re- 
move background, the invariant mass of the good photon 
and the missing momentum must be consistent with that 
of the TT*^, and the invariant mass of the charged tracks 
from the J/tp decay must be consistent with that of the 
p. The photon detection efficiency is then the fraction 
of actual photons matched in direction to the predicted 
photon. On average, the efficiency difference between 
data and Monte Carlo simulation is less than 1%. 

The missing 7r° method uses a sample of V' ^ 
7r°7r° J/V", J/i/j events, in a similar way to the 

first method. Events with two charged tracks and at 



least two photons are required. The invariant mass of 
the charged tracks must be consistent with the mass of 
the J/ip, and the sum of the momenta of the two pho- 
tons must be greater than 300 MeV/c. Since the two 
pions are anti-correlated, by requiring a tt'^ with momen- 
tum larger than 300 McV/c and using energy momentum 
conservation, the energy and direction of the remaining 
soft 7r° can be predicted. To ensure a clean sample, the 
invariant mass of the two photons of the selected pion 
and the mass recoiling against the J/ip and the selected 
7r° must be consistent with that of the 7r°. The number 
of reconstructed tt'^s, which match in direction and have 
two photon invariant mass consistent with the tt*^ mass, 
yields the 7r° detection efficiency. The difference in ef- 
ficiency between the data and Monte Carlo is less than 
(1.5 ± 0.5)%; hence the uncertainty of photon detection 
is less than 1%. 

The third method, the 7r° decay angle method [l^ . 
utilizes a J/-0 — ?> /9°7r° sample. Since the 7r° is from a 
two-body decay, the momentum of the tt^ is known, and 
the energy of the lower energy daughter photon is given 

by 2Eiow = \J P^o + M^o — PttO cos 9, where 9 is the angle 

of the photon in the tt" rest frame with respect to the di- 
rection of the TT*^ in the J/ ip rest frame. For spin-zero par- 
ticles like the 7r° , the cos 6 distribution is flat in the range 
of [0,1]. Hence the Eiow energy spectrum is intrinsically 
fiat, and the deviation from fiatness measures the pho- 
ton inefficiency as a function of energy. The distribution 
shows that the efficiency is low at low energy but plateaus 
starting from 0.1 GcV. A comparison of the photon en- 
ergy spectrum shape between data and Monte Carlo sim- 
ulation shows that the average difference is 0.6%. How- 
ever, this measurement only provides a relative efficiency 
as a function of photon energy. The absolute efficiency 
in the plateau region above 0.1 GeV can be determined 
in data and simulation using high energy electrons from 
radiative Bhabha events. The electromagnetic showers 
of electrons and photons in the crystal calorimeter are 
the same for E > 0.2 GeV. The detection efficiency for 
electrons entering the EMC with E > 0.2 GeV is essen- 
tially 100% in both data and Monte Carlo simulation, 
which indicates no significant systematic uncertainty in 
the simulated efficiency of the EMC. 

The above three methods may be affected by pho- 
ton conversions, mainly in the Beryllium beam pipe 
and inner part of the MDC. A study using samples of 
e+e~ — >■ 77,7 — >■ e~^e~ events in which the converted 
photons are explicitly reconstructed indicates systematic 
efficiency differences due to material in the interaction 
region between data and simulation are negligible. 

Although each of the above three methods suffers from 
different shortcomings, such as the resolution and track- 
ing efficiency dependence, they all give consistent results 
within 1%, showing that the photon efficiency uncer- 
tainty is less than 1%. 
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B. Tracking efficiency 



TABLE II: Sources of systematic errors (%). 



The tracking efBcicncics for soft and hard pions arc 
studied with ti^-k^ J J /4> —J- — e,ii) 

and J/ip — > pTT — >■ Tr+Tr^TT*^ event samples, respectively. 
The transverse momentum for a soft pion is less than 
400 MeV/c. The tracking efficiency is calculated with 
e = N full /Nail, where Nfuii indicates the number of 
events of 7r+7r~Z+Z~(7r+7r~7r°) with all final tracks re- 
constructed successfully; Naii indicates the number of 
events with one or both charged pion tracks success- 
fully reconstructed in addition to the lepton pair (tt") for 

V^' TT+TT" J/?/), J/-0 I'^l" {Jj'^ -> pTT -> 7r+7r~7r*'). 

In addition, we require that the direction of the miss- 
ing momentum should be within the MDC coverage. 
The missing momentum is calculated using the recon- 
structed Icpton-pair (tt") and one of the reconstructed 
pions for the -0' — >■ tt+tt" J/^, J/?/) — >■ (•//V' — >■ Ptt ^• 
7r+7r~7r°). A very clean soft pion sample is selected by 
doing a kinematic fit with the lepton pair constrained to 
the J/?/' mass, and background from J/^ — >■ tt+tt" is re- 
jected using muon counter information and is negligible. 
For the hard pion selected from J /tp 7r+7r~7r'', the pu- 
rity of the sample is more that 98%; the small background 
is from J/-0 — ?> K*K — >■ K^K~-k^ due to misidcntifica- 
tion of the kaon as a pion for the reconstructed track. 
The differences for the soft and hard pion tracking effi- 
ciencies between the data and MC are both estimated to 
be 2%, which is taken as the pion tracking uncertainty. 

The kaon tracking efficiency is determined with a sam- 
ple of J/V- K*{%mfKl + c.c ^ KlK+TT- + c.c ^ 
i^r+TT^TT+TT^ -f i4r~7r+7r+7r~ events. The tracking effi- 
ciency is calculated in the same way as the pion track 
efficiency. Here Nfuii is the number of events with a 
matched kaon track in addition to the three tracks iden- 
tified as pions, and Naii is the number of events with or 
without a matched kaon and three tracks identified as 
pions. A clean sample is selected by using a second ver- 
tex fit to Kg 1 and setting stringent mass windows for the 
Ks and iir*(892)°. The difference in the kaon tracking 
efficiency is about 2% between the data and MC, which 
is taken as the uncertainty of kaon tracking efficiency. 



C. Efficiency for particle identification 

The efficiencies for pion and kaon PID are obtained 
with J/ip 7r+7r~7r° and K^K^-k^ control samples, 
respectively. Samples with backgrounds less than 1% 
are selected by using a narrow tt" mass window, and 
requiring one track be identified as a pion (kaon) for 
J/ip 7r+7r-7r° {K+K'tt^) based on the TOF and 
dE/dx information. The PID efficiency is calculated with 
e{PID) = No/{N' + No), where Nq {N') denotes the 
events with the other track identified (not identified) as 
a pion or kaon. The differences between data and MC 
for the pion and kaon PID efficiencies arc about 2%, and 
2% is taken as the systematic error. 
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2.1, 2.0, 2.5 


1.0, 0.9, 1.5 


1.0, 0.5, 1.0 


Signal shape 


0.8 


0.8 


0.8 


No. of ijj' evts 


3.8 


3.8 


3.8 


^ 7Xcj) [15j 


3.3, 4.4, 4.0 


3.3, 4.4, 4.0 


3.3, 4.4, 4.0 


V decay [15] 


1.0 




0.8 


Total 


8.8, 8.8, 9.3 


9.5, 9.7, 10.0 


9.3, 9.4, 9.8 



D. Uncertainties on the selection efficiencies 



The uncertainties on the selection efficiencies listed in 
Table|n]are considered as the systematic errors due to the 
selection criteria. Control samples such as J/i/' — > p°7r°, 
J/V' — >■ ipVj and J/ip LOT] are used to determine the 
efficiency difference between data and MC for each se- 
lection criteria. For ip' — >■ j"f4>, the uncertainty for the 
?7 veto is estimated to be 1.9%, and that for the selec- 
tion of the (j) signal is 0.5%. In -0' jjp, the dominant 
sources are due to the electron and muon track rejection, 



the cos requirement, the tt /ry veto, and the selection 
of the p° signal, which are 4.5%, 0.9%, 1.2%, and 1.4%, 
respectively. For this channel, the uncertainty due to the 
muon track veto for both the 7r+ and tt" (4.5%) is deter- 
mined from a study of a very pure J/ip — > p^n'^ sample. 
In 4'' ~> 77^ J the uncertainty due to the background 
rejection and w signal selection is 1.4%. The total uncer- 
tainties due to selection criteria are 2.0%, 5.0%, and 1.4% 
for xjy — > 770, ip' ^ 77P°7 and i/j' — > 77a;, respectively. 

Systematic errors in the fit to the ^hV mass distribu- 
tion originate from the uncertainties in the parameteri- 
zations for the signal and background shapes. We obtain 
the uncertainties due to the background shape by chang- 
ing the order of the polynomial function from 2 to 3 in the 
fit. and those for the signal shapes arc estimated by con- 
volving a smearing Gaussian function to account for the 
difference in the resolution between data and MC sim- 
ulation. Systematic errors associated with the bin size, 
fitting range, and vector meson side-band regions are es- 
timated by changing the bin size, the fitting range, and 
the side-band regions. 

All errors are summarized in Table|lTl Finally, the total 
systematic error varies from 8.8% - 10.0% depending on 
the final state as summarized in Table |TT1 
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V. HELICITY AMPLITUDE ANALYSIS 

In Xci — >■ decays, the final state is a superposition 
of longitudinal (A — 0) and transverse (A = ±1) polar- 
izations. The angular distribution is 

— ^ cx (1 - /t) cos2 e + i/T sin^ e, 
r d cos 2 

where /t = |AtP/(|AtP + |^lP) is the transverse polar- 
ization fraction in the decay and Al and At are the lon- 
gitudinal and transverse polarization amplitudes, respec- 
tively, and Q is the angle between the vector meson flight 
direction in the Xci rest frame and either the /K^ di- 
rection in the p'^ /(f> rest frame or the normal to the w 
decay plane in the uj rest frame. By performing a like- 
lihood fit to the angular distributions of the vector me- 
son decays, we can determine the transverse polarization 
fraction Jt- We account for the different reconstruction 
efficiencies for the longitudinally and transversely polar- 
ized events using fr = f^^y{R+{l-R)f^^''), where f^^"" 
is the fraction of signal from transversely polarized signal 
events in data, and R is the ratio of the longitudinal and 
transverse signal efficiencies. 

In the likelihood fit, we take events in the Xci signal 
region which is defined as 3.49 < M^^y < 3.52 GeV/c^. 
The signal shapes for the longitudinal and transverse 
component are obtained from MC simulations. For the 
background shapes, two sources are considered: one is the 
vector-meson sideband background, which is normalized 
according to the number of sideband events in the Xci 
signal region; the other background considered is from 
MC simulated inclusive tp' decay events (not including 
signal events) that satisfy the selection criteria and have 
an invariant mass in the Xci signal region, which are nor- 
malized according the number of polynomial events used 
in the Xci mass fit. The cos Q distributions are fitted with 
the combined backgrounds and MC simulated transverse 
and longitudinal signal shapes. The total signal yield 
and transverse polarization fraction are floated in the fit. 
The fitted results are shown in Fig. [3] The values of the 
fraction of the transverse component are 0.291q Jj^"'^*' 



CLEO-c's A±/Ao have been used to determine fx as: 

072+° °'*i+° °°2 for y 1 TO and 32+o i7+o.05 r 

iJ-" ' ^-0.031-0.019 ^ Acl ^ ;Pi aim u.o^_q -^j^_g j^-^ loi 

Xci 7^^- The results are consistent within 2a. 



for Xci ^ 70, 0.158 ± 0.034t° °M for Xa 

7-1-0.090+0.044 



IP 



0.09 

and 



0.247_g;o87_o.o26 fo^ Xa ^ 7^, where the first errors arc 
statistical from the fit and the second ones arc the sys- 
tematic errors. 

Since fx is a ratio, many systematic errors cancel out, 
and only the effects due to binning of the cos Q distribu- 
tions and the parameterization of the background shape 
(estimated by assuming the backgrounds (except for the 
vector-meson sideband background) contribute entirely 
to either the longitudinal or the transverse component) 
are considered here. The uncertainties of the binning of 
cos 9 are estimated to be ^ig ^^, -ae*^, and 11q[1%, and 
the parametrization of the background, which is the dom- 
inant systematic error, is estimated to be ^9^2'^ 

, and tio:5% for Xa 70, Xa ~^ 1P°- and Xci ^ joj, 
respectively. In order to compare with CLEO-c's results. 




■0.4 -0.2 „ 0.2 0.4 0.6 0.8 
COS0 



FIG. 3: cos distributions and the fit for (a) xa — >■ 70, (b) 
Xci — !■ 7p", and (c) xa — Ji^- Dots with error bars are data; 
the solid histogram is the fit; and dashed histogram is the sum 
of the sideband background and the polynomial background. 



VI. FINAL RESULTS AND DISCUSSION 

In summary, we present the measurements of radiative 
decays of Xcj to light vector mesons. We find B{xci 
-ip^) = (228±13±16) X 10-6 and B{xci ^ 7^) = (69.7± 
7.2 ± 5.6) X IQ-^, which agree with the results from the 
CLEO experiment Q . We observe Xci l4> for the first 
time, and find B{xci ^ 70) = (25.8 ± 5.2 ± 2.0) x 10"^. 
Upper limits at the 90% confidence level on the branch- 
ing fractions for Xco and Xc2 decays into these final states 
arc determined. The final results are listed in Table IIIII 
The theoretical predictions for B(xci iV) including 
the hadronic loop contribution in pQCD calculation Q 
are consistent with our measurements within errors. In 
addition, the fraction of the transverse polarization com- 
ponent of the vector meson in Xci — ^ iV decay is stud- 
ied. Our measurements of the polarization of the vector 
mesons indicate that the longitudinal component is dom- 
inant in Xci — ^ iV decay, as expected for an axial-vector 
particle radiative decaying into a vector ((/>, p^ , and u]) in 
the framework of the vector dominance model taking into 
account the Landau- Yang theorem |2(|,i2lj . This observa- 
tion may aid future development on the QCD calculation 
of the partial waves in the Xci ~^ iV decay. 
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TABLE III: Results on Xcj ~^ "fV . The upper limits are set 
at the 90% C.L. 



Decay 




No. of 


Eff. 


Syst. 


Br. 


Stat. 


mode 




evts. 


(%) 


err.(%) 


(io-«) 


sign. 


Xco 




15.0 ±6.6 


32.4 


8.8 


<16.2 




Xcl 


70 


42.6 ±8.6 


34.6 


8.8 


25.8 ±5.2 ±2.3 


6cr 


Xc2 


70 


4.6 ±4.9 


32.6 


9.3 


<8.1 




XcO 


7P» 


6± 12 


22.6 


8.1 


<10.5 




Xcl 


7P" 


432 ± 25 


19.4 


8.3 


228 ± 13 ± 22 


> lOcr 


Xc2 


7P° 


13 ± 11 


15.7 


8.7 


<20.8 




Xco 


7CJ 


5± 11 


18.6 


9.3 


<12.9 






70; 


136 ± 14 


22.7 


9.4 


69.7 ±7.2 ±6.6 


> 10(7 


Xc2 


7tj 


1±6 


19.2 


9.8 


<6.1 
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